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A  nonlinear,  compressible,  spectral  collocation  code  has  been  developed  to  exam in 
gravity  wave  breaking  and  instability  processes  in  two  and  three  spatial  dimen_i.; 
Initial  studies  have  demonstrated  that  the  preferred  mode  of  instability  within  a 
high-frequency  gravity  wave  is  a  convective  instability  comprised  of  counter-rota 
ing  vortices  aligned  transverse  to  the  direction  of  wave  propagation  (a  horizonta 
wavenumber  normal  to  that  of  the  gravity  wave).  Thus,  wave  instability  is  inhere 
three-dimensional,  and  two-dimensional  models  are  unlikely  to  adequately  describe 
either  the  physics  of  wave  breaking  or  the  implications  for  wave  transports  and 
eddy  mixing.  A  parallel  effort  has  emphasized  the  statistical  effects  of  wave 
interactions  and  dissipation  processes  and  developed  a  new  spectral  parameter  iz.ft 
of  gravity  wave  transports  of  energy  and  momentum  and  their  atmospheric  effects. 
This  scheme  relies  on  the  approximately  universal  spectral  shape  of  the  gravity- 
motion  field  throughout  the  atmosphere  to  assess  the  potential  for  wave  transport, 
and  variations  with  background  wind  and  stability  profiles. 
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2.  Research  Objectives 

The  goals  of  this  research  were  to  address  two  sets  of  problem*  relating  to  grav¬ 
ity  wave  effects  in  the  atmosphere.  One  approach  was  intended  to  apply  analytic 
and  pseudo-spectral  modeling  techniques  to  studies  of  wave  excitation,  instability, 
dissipation,  and  interaction  processes  thought  to  play  important  roles  in  the  cir¬ 
culation.  structure,  and  variability  of  the  atmosphere  at  lower  and  upper  levels. 
A  second  effort  was  intended  to  apply  our  current  understanding  of  the  spectral 
character  of  the  atmospheric  motion  field  to  the  development  of  a  parameterization 
scheme  for  the  most  important  wave  effects  in  large-scale  models  of  the  atmosphere. 
Substantial  progress  has  been  made  in  both  areas  and  is  discussed  further  below. 

3.  Research  Progress  to  Date 

A  nonlinear,  compressible,  spectral  collocation  code  has  been  developed  to 
examine  gravity  wave  breaking  and  instability  processes  in  two  and  three  spatial 
dimensions.  Initial  studies  have  demonstrated  that  the  preferred  mode  of  instabil¬ 
ity  within  a  high-frequency  gravity  wave  is  a  convective  instability  comprised  of 
counter-rotating  vortices  aligned  transverse  to  the  direction  of  wave  propagation  (a 
horizontal  wavenumber  normal  to  that  of  the  gravity  wave).  Thus,  wave  instability 
is  inherently  three-dimensional,  and  two-dimensional  models  are  unlikely  to  ade¬ 
quately  describe  either  the  physics  of  wave  breaking  or  the  implications  for  wave- 
transports  and  eddy  mixing.  These  studies  are  described  in  .*  -cries  of  three  papers 
that  have  been  submitted  to  J.  Geophys.  Res.  Additional  numerical  studies  have 
addressed  the  effects  of  a  transverse  shear  on  wave  instability  and  are  now  being 
applied  to  shear  instability  problems  of  relevance  to  mixing  and  transports  in  the 
atmosphere. 

A  parallel  effort  has  emphasized  the  statistical  effects  of  wav.-  interactions  and 
dissipation  processes  and  developed  a  new  spectral  parameterization  of  gravity  wave 
transports  of  energy  and  momentum  and  their  atmospheric  effects.  This  scheme 
relies  on  the  approximately  universal  spectra1  shape  of  the  giavity  wave  motion 
field  throughout  the  atmosphere  to  assess  the  potential  for  wave  transports  and 
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variations  with  background  wind  and  stability  profiles,  Tm  is  presently 

working  in  the  NCAR  TGCM,  is  being  tested  by  a  number  of  other  modeling 
groups,  and  is  described  in  a  series  of  three  papers  presently  in  press  in  J.  Atmos. 
Sci. 

Finally,  we  have  performed  a  series  of  studies  of  specific  ginvity  wave  sources 
or  source  processes  using  both  theoretical  and  data  analysis  techniques.  The  latter 
have  shown  the  dominant  sources  of  small-scale  wave  motions  to  be  orogranhy, 
convection,  and  wind  shear.  Analytic  studies  have  examined  the  excitation  of 
gravity  waves  by  geostrophic  adjustment  of  the  jet  stream  and  by  eclipse  forcing. 
Those  papers  citing  AFOSR  support  are  listed  following  this  report. 

4.  Future  Research  Efforts 

Our  future  efforts  under  AFOSR  support  will  continue  to  address  the  physics 
and  consequences  of  wave  propagation  and  instability,  with  emphasis  on  both  large- 
and  small-scale  effects.  At  small  scales,  our  simulations  will  address  the  transition 
from  laminar  to  turbulent  flow,  the  intermittency  and  evolution  of  turbulence  un¬ 
der  the  influences  of  stratification,  and  the  induced  transports  of  momentum.  heat, 
and  constituents.  Such  issues  are  important  in  understanding  the  effects  of  wave 
instability  on  the  thermal  and  constituent  structures  of  the  lower  and  middle  at¬ 
mosphere.  At  larger  scales,  our  efforts  will  address  the  constraints  imposed  by 
instability  processes  on  wave  amplitudes  and  transports,  the  spectral  evolution 
of  the  motion  field,  the  effects  of  wave- wave  and  wave-mean  flow  interactions  on 
wave  propagation,  and  wave  influences  on  the  large-scale  circulation  and  thermal 
structure  via  their  fluxes  of  energy  and  momentum. 

Specific  efforts  will  include  1)  studies  of  wave  breaking  for  lower-frequency 
motions  which  may  favor  other  modes  of  instability,  2)  simulations  addressing 
instability  processes,  effects,  and  variability  in  a  multiple- wave  environment,  3) 
studies  of  the  effects  of  instability  and  wave-wave  interactions  on  the  spectral  evo¬ 
lution  of  gravity  waves  in  the  atmosphere,  4)  and  wave  and  instability  responses 
to  specific  sources  or  propagation  environments,  such  as  mountain  wave  instability 
and  transience  and  the  influences  on  a  gravity  wave  spectrum  of  a  tidally-varying 
background  wind  and  thermal  field. 

These  studies  will  take  advantage  of  new  versions  of  our  spectral  collocation 
code  that  will  employ  the  anelastic  approximation  in  order  to  peiinit  higher  resolu¬ 
tion  and  longer  integration  times,  relative  to  our  previous  compressible  simulations. 
Present  models  using  Fourier  basis  functions  in  the  horizontal  and  Chebvshev  poly¬ 
nomials  in  the  vertical  will  be  applied  to  studies  of  shear  instability,  multiple  wave 
saturation  and  spectral  evolution,  and  turbulence  production  and  decay.  Like  grav¬ 
ity  wave  breaking,  shear  instabilities  are  very  important  to  our  understanding  ot 
transports  and  mixing  in  the  atmosphere  because  the  conditions  leading  to  insta¬ 
bility  are  quite  different  than  those  within  a  breaking  gravity  wave.  As  such,  we 
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anticipate  that  their  implications  for  turbulent  transports  may  also  be  very  differ¬ 
ent  than  wave  breaking.  Yet  there  have  been  no  studies  to  date  that  have  addressed 
the  full  three-dimensional  evolution  of  such  instabilities.  High-resolution  numerical 
models  are  also  likely  the  only  means  of  studying  the  full  nonlinear  evolution  of  a 
wide  spectrum  of  wave  motions  and  of  testing  the  various  theories  that  have  been 
advanced  to  try  and  account  for  observed  spectral  features  of  tin-  atmospheric  wave 
field.  Computation  of  these  wave  and  instability  evolutions  to  late  times  will  also 
address  the  morphology  of  turbulence  in  stratified  flows  and.  address  problems  of 
fundamental  importance  for  tracer  transports  and  measurement  techniques  relying 
on  turbulent  backscatter. 

Another  version  of  our  spectral  collocation  code  using  Cheby-hev  polynomials 
in  the  horizontal  and  vertical  will  be  applied  to  problems  such  as  mountain  wave 
structure  and  instability  for  which  inflow  and  outflow  condition-  may  nor  be  pe¬ 
riodic.  These  flows  depend  strongly  on  instability  processes  and  variability,  and 
accurate  descriptions  of  the  nonlinear  interactions  involved  are  essential  to  under¬ 
standing  the  process  as  a  whole  and  its  implications  for  the  atmosphere  at  larger 
scales. 

Our  parameterization  efforts  will  aim  to  further  refine  the  gravity  wave  param¬ 
eterization  scheme  developed  under  current  AFOSR  support  and  now  being  tested 
and/or  implemented  in  a  number  of  mechanistic  and  GCM  models  by  rations 
groups.  This  parameterization  scheme  was  designed  with  the  intent  of  describing 
statistically  the  responses  of  the  gravity  wave  spectrum  to  varying  environme:/ 
and  its  influences  on  the  local  mean  flow  that  result  from  wave  saturation  tm  i 
filtering.  As  such,  it  responds  in  an  averaged  manner  to  variations  in  mean  wind- 
that  is  in  approximate  agreement  with  the  limited  observations  that  are  available 
The  scheme  also  has  the  ability  to  describe  the  local  responses  to  local  gravity  wave 
sources  posed  at  discrete  locations  and  times.  However,  the  scheme  does  not .  at 
present,  have  a  procedure  for  distributing  these  local  source  effects  geographically 
with  increasing  height.  At  greater  altitudes,  it  also  is  important  to  include  a  nman- 
of  describing  the  anticipated  temporal  variations  in  wave  dissipation  and  ir.du*  •  1 
diffusion  in  order  to  account  for  the  tidal  modulation  of  these  pi ocesses  and  tie  .: 
influences  in  turn  on  the  transports  of  heat  and  trace  constituents.  Additional  < 
forts  will  advance  the  present  scheme,  test  its  responses  to  variou-  ■  voonmen's 
1-D  and  2-D  mechanistic  models  as  well  as  in  various  GCM's.  su>  1.  the  TGCM 
at  NCAR,  which  is  now  using  the  present  version  of  the  param-  t<  ;  n  ation, 
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